Effect of spin-conserving scattering on Gilbert damping in ferromagnetic 

semiconductors 
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The Gilbert damping in ferromagnetic semiconductors is theoretically investigated based on the 
s-d model. In contrast to the situation in metals, all the spin-conserving scattering in ferromagnetic 
semiconductors supplies an additional spin relaxation channel due to the momentum dependent 
effective magnetic field of the spin-orbit coupling, thereby modifies the Gilbert damping. In the 
presence of a pure spin current, we predict a new contribution due to the interplay of the anisotropic 
spin-orbit coupling and a pure spin current. 
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The ferromagnetic systems have attracted much at- 
tention both for the abundant fundamental physics and 
promising applications in the past decade. ^'^ The study 
on the collective magnetization dynamics in such sys- 
tems has been an active field with the aim to control 
the magnetization. In the literature, the magnetization 
dynamics is usually described by the phenomenological 
Landau-Lifshitz-Gilbert (LLG) equation)^ 



eff 



(1) 



with n denoting the direction of the magnetization. The 
first and second terms on the right hand side of the equa- 
tion represent the precession and relaxation of the mag- 
netization under the effective magnetic field Hog , respec- 
tively. The relaxation term is conventionally named as 
the Gilbert damping term with the damping coefficient 
a. The time scale of the magnetization relaxation then 
can be estimated by l/(a7iJeff)j^ which is an important 
parameter for dynamic manipulations. The coefficient a 
is essential in determining the efficiency of the current- 
induced magnetization swiching, and experimental deter- 
mination of a has been carried out intensively in metals^ 
and magnetic semiconductors.- 

To date, many efforts have been made to clarify the 
microscopic origin of the Gilbert dampingi^r— Kohno 
et al,^ employed the standard diagrammatic pertur- 
bation approach to calculate the spin torque in the 
small-amplitude magnetization dynamics and obtained a 
Gilbert torque with the damping coefficient inversely pro- 
portional to the electron spin lifetime. They showed that 
the electron-non-magnetic impurity scattering, a spin- 
conserving process, does not affect the Gilbert damping. 
Later, they extended the theory into the finite-amplitude 
dynamics by introducing an SU(2) gauge field2 and ob- 
tained a Gilbert torque identical to that in the case of 
small-amplitude dynamics.^ In those calculations, the 
electron-phonon and electron-electron scatterings were 
discarded. One may infer that both of them should be 
irrelevant to the Gilbert damping in ferromagnetic met- 
als, since they are independent of the electron spin re- 



laxation somewhat like the electron-non-magnetic impu- 
rity scattering. However, the situation is quite different 
in ferromagnetic semiconductors, where the spin-orbit 
coupling (SOC) due to the bulk inversion asymmetrjii^ 
and/or the structure inversion asymmetry^** presents a 
momentum-dependent effective magnetic field (inhomo- 
geneous broadeningi^) . As a result, any spin-conserving 
scattering, including the electron-electron Coulomb scat- 
tering, can result in a spin relaxation channel to affect the 
Gilbert damping. In this case, many-body effects on the 
Gilbert damping due to the electron-electron Coulomb 
scattering should be expected. Sinova et ali^ studied the 
Gilbert damping in GaMnAs ferromagnetic semiconduc- 
tors by including the SOC to the energy band structure. 
In that work, the dynamics of the carrier spin coherence 
was missed.— The issue of the present work is to study 
the Gilbert damping in a coherent frame. 

In this Report, we apply the gauge field approach to 
investigate the Gilbert damping in ferromagnetic semi- 
conductors. In our frame, all the relevant scattering pro- 
cesses, even the electron-electron scattering which gives 
rise to many-body effects, can be included. The goal 
of this work is to illustrate the role of the SOC and 
spin-conserving scattering on Gilbert damping. We show 
that the spin-conserving scattering can affect the Gilbert 
damping due to the contribution on spin relaxation pro- 
cess. We also discuss the case with a pure spin current, 
from which we predict a new Gilbert torque due to the 
interplay of the SOC and the spin current. 

Our calculation is based on the s-d model with itiner- 
ant s and localized d electrons. The collective magnetiza- 
tion arising from the d electrons is denoted by M = M^n. 
The exchange interaction between itinerant and local- 
ized electrons can be written as Hgd = M J dr{n ■ cr), 
where the Pauli matrices cr are spin operators of the 
itinerant electrons and M is the coupling constant. In 
order to treat the magnetization dynamics with an ar- 
bitrary amplitude," we define the temporal spinor oper- 
ators of the itinerant electrons a{t) — {af{t),ai{t))'^ in 
the rotation coordinate system with ^ (i) labeling the 
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spin orientation parallel (antiparallel) to n. With a uni- 
tary transformation matrix U(t), one can connect the 
operators a-t'(4.) with those defined in the lattice coor- 
dinate system c^(i) by a{t) — U{t)c. Then, an SU(2) 
gauge field = ~iU{t)\dfJJ{t)) = A^{t) ■ a should 

be introduced into the rotation framework to guarantee 
the invariance of the total Lagrangian.-^ In the slow and 
smooth precession limit, the gauge field can be treated 
perturbativelyi^ Besides, one needs a time-dependent 
3x3 orthogonal rotation matrix Jl{t), which obeys 
U^(tU — TZcr, to transform any vector between the two 
coordinate systems. More details can be found in Ref. 
0. In the following, we restrict our derivation to a spa- 
tially homogeneous system, to obtain the Gilbert damp- 
ing torque. 

Up to the first order, the interaction Hamiltonian due 
to the gauge field is Ha ~ X^k -^o • aj^crak and the spin- 
orbit couping reads 



c' crc 



k • al^cru]^, 



(2) 



where pk represent the itinerant electron density matrices 
defined in the rotation coordinate system. The coherent 
term can be written as 



^tPi^Lh = -i[A- 



^HF, PkJ 



(4) 



Here [, ] is the commutator and A{t) = Ao(t) + Mz with 
Aq and Mz representing the gauge field and effective 
magnetic filed due to s-d exchange interaction, respec- 
tively. Shf is the Coulomb Hartree-Fock term of the 
electron-electron interaction. 9tPk ^ and dtOh-l , in 
Eq. ^ include all the relevant spin-conserving and spin- 
fiip scattering processes, respectively. 

The spin-flip term dtpk\^^ results in the damping ef- 
fect was studied in Ref. ^ Let us confirm this by 
considering the case of the magnetic disorder = 

Us J2j ' 0'''o-ciS{r — Rj). The spin-fiip part then reads 



with h = T^h. Here, we take the Planck constant H = 1. 
We start from the fully microscopic kinetic spin Bloch 
equations of the itinerant electrons derived from the non- 
equilibrium Green's function approach , ^^'^^ 



dtPk = dtPk 



coh 



dtP\^\ 



dtPk 



f 



Q l/ O |/(0) 



with dtPkl standing for the i-th order term with re- 



dtPk 



/(I) 



(5) 



(3) spect to the gauge field, i.e.. 



dtP^Cl = E 'T"p>(t)r,,a"T„,p<(t)5(6k,,, -ek,.)-(>o<)+H.c., (6) 

kl77iJ)2 

dtPk\i^, = ^^e-Z^-r^^W E '^"<Wr..^''r,,p<W- <S(ek,,, -ek,.)-(>o<)+H.c., (7) 

ki?)i?72 



where T,j{i,j) = S^^5^^ for the spin band rj. Here 
p^ = 1 — Pk, Pk ~ Pk- ( is obtained by inter- 
changing > and < from the first term in each equation, 
e'-''^ is the Levi-Civita permutation symbol. The gauge 
field term, dtPk\^^^^ , results from the spin correlation of 
a single magnetic impurity at different times It induces 
a spin polarization proportional to z x (t) which gives 
a Gilbert torque. The damping coefficient is inversely 
proportional to the spin relaxation time detemined by 

the spin-flip scattering i9tPk|gj|^j • The spin-fiip scattering 
term in Eq. ^ thus reproduces the result of Ref. |9|. 

We now demonstrate that the Gilbert damping torque 
arises also from the spin-conserving scattering. For 
the discussion of the spin-conserving term, it is suffi- 

I / 

cient to approximate the spin-flip term as dtpk\^^^^, — 
— (Pk — Pk )/■''« J with p^ representing the instantaneous 
equilibrium distribution (i.e., p^ is pk without the gauge 



field and P^)- Equation ([3]) then reads 

9tPk = -«[-4-cr,pk] -^[^hk •cr,pk] 

+ 9tPk|lt-(Pk-p^)/T. + PL (8) 

Here, we add an additional term, P^, to describe the 
source of a pure spin current due to the magnetization 
dynamic pumping^ or electrically injectioni^i^ in order 
to discuss the system with a pure spin current. We ne- 
glect the Coulomb Hartree-Fock effective magnetic field 
since it is approximately parallel to the s-d exchange field, 
but with a smaller magnitude. 

By averaging density matrices over the momentum 
direction, one obtains the isotropic component pi k — 
J ^§^Pk- The anisotropic component is then expressed 
as Pa.k = Pk ~ Pi.k- It is obvious that this anisotropic 
component does not give any spin torque in the absence 
of the SOC, since J^k^^i^Pa-k) = 0. Below, it is shown 
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that this component leads to the damping when coupled 
to the spin-orbit interaction. 

By denoting the isotropic component of the equi- 
librium part (p^) as pl^ and representing the non- 
equilibrium isotropic part by 5pi_k = Pi.k — Plk^ 
write the kinetic spin Bloch equations of the non- 
equilibrium isotropic density matrices Spi^k and those of 
the anisotropic components pa;k as 



dtpi,k 

dtPa,k 



i[A- cr,Spi^k] 
-i[Ao-(T,pl^], (9) 

dtPa,k\l^^^ - <T,Pa,k] -j[5hk •cr,'5pi,fe] 



i[ihk ■ (T, pa,k] + i[hhk ■ cr, pa^] + Pk(lO) 



respectively. The overline in these equations presents a 
angular average over the momentum space. 

We further define p^^^ as the anisotropic density in the 
absence of the gauge field, Aq. As easily seen, it vanishes 
when — 0. The anisotropic component involving the 

gauge field is denoted by p^^]^ = Pa.k — pi^k- Equation 
(1101) is expressed in terms of these components as 



^1 (0) 



a (1) 



-^[M■cT,p'"l]+^tp'[ 



k I scat 



5hk-cr,pi°^]+z[ihk-cr,p^'; 



(0)i 



kJ' 



(11) 



^[^h^■cT,pi'h-^[Ao■a,pl'i] 



(0)i 



(12) 



Within the elastic scattering approximation, the 
electron-phonon scattering as well as the electron-non- 
magnetic impurity scattering can be simply written as 

J2i,mP''a,k,im^'-"^/^'^ whcrc the density matrices are ex- 
panded by the spherical harmonics functions 1^™, i.e., 

Pal.im = I ^Pa lYim- n is the effective momentum 
relaxation time. The exact calculation of the Coulomb 
scattering is more complicated. Nevertheless, one can 
still express this term in the form of pi^kZ-^klp), where 
^k is a function of the density matrices^! and reflects 
many-body effects. For simplification, we just introduce 
a uniform momentum relaxation time r* in the following 
calculation. Expanding Eq. by the spherical har- 
monics functions, one obtains 

^tPilim = 'i\-^-<^^P^al,lm]-^\\^k.l^-(r,5p,,k\ 



(1) 



^^,13) 



where the expanssion coefficient of any term /k is ex- 
pressed as fk.im = / ^IkYim- In the strong scattering 
regime, i.e., 4r >• M and \ » Ihkl, the first and fourth 
terms are much smaller than the last term, hence can be 
discarded from the right side. By taking the fact that 
the time derivative is a higher order term into account, 
one also neglects dtp^^\ ;,„. The solution of Eq. ([TB]) can 



be written as 

pi!L™ = -*r;{[ihfc,,,„.^,5p,,,] + [Ao-^,pi°l.^J}. (14) 

Substituting it into Eq. (jl4p and rewriting the equation 
in the leading order, one obtains 



^tp^,k = -i[A ■ (T, 6p,^k] - ^[hk • cr, p^^]^] - i[Aa ■ cr, pl^] 
-4- ^k,im ■ cr, [hfc,/m • cr, dpi^k 



(15) 



The third term on the right hand side of the equation is 
proportional to the second order term of the SOC, which 
gives the spin dephasing channel due to the D'yakonov- 
Perel' (DP) mechanism .=22, This term can be expressed 
by TjjpSpi^k with T^p standing for the spin dephas- 
ing rate tensor, which can be written as {T^p)ij — 
J2i,,n{Ti'ii^k,ini)^Sij-hl^i^hi^i^)) by performing the en- 
semble averaging over the electron distribution. In the 
following, we treat tdp as a scalar for simplification and 
the total spin lifetime is hence given by 



= 1/0 



DP 



(16) 



Similar to the previous procedure, we discard dtPi^k in 
Eq. and obtain 



i[A-cr,Spi,k]+Spt,k/Tr = -i[^hk • (T,pl l.]-i[Ao-cr, p'l ^.]. 

(17) 

By taking Ss, - i Ek Tr(<T^P,,/c), = hT.k'^^i'^Plk) 



and s'-°l 



iTr(crp[^'^^), one can write the solution as 



V + 2TrA X V + (v . ^) ^ ^ 



l + 4|^|2r2 



(18) 



Tr J2k ^k X s^^l^. sf is just the equilibrium 



where v — s. 

spin density , which is parallel to the magnetization, i.e., 
s,^ = sfz. Now, we pick up the transverse component in 
the form of z x Ag , Ss-^ , since only this component results 
in a Gilbert torque of the magnetization as mentioned 
above. We come to 



(19) 



with Tox 

;oc 
obtains 



1/(2M). By transforming it back to the lat- 

^dtn^ one 



tice coordinate system with TZ{z x Aq ) = ^dfU — 



-V^{dtn)T^^Tr/{T^ 



(20) 



This nonequilibrium spin polarization results in a spin 
torque performed on the magnetization according to T = 
— 2Mn X Ss, i.e.. 



■0^(11 X dtn)TcycTr / [tI^ + r^x)- 



(21) 



Compared with Eq. ([1]), the modification of the Gilbert 
damping coefficient from this torque is 



(22) 
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We fisrt discuss the case without the source term of 
the spin current. In this case, the anisotropic component 
Pa°k vanishes and Vz = s^. We see that the Gilbert damp- 
ing then arises from I/t^ [Eq. ([T6|]. i.e., from both the 
spin-flip scattering and the DP mechanism/^ Our main 
message is that this DP contribution is affected by the 
spin-conserving scattering processes such as the electron- 
electron interaction and phonons. The temperature de- 
pendence of the Gilbert damping and the current- induced 
magnetization switching can thus be discussed quantita- 
tively by evaluating r^. We note that our result reduces 
to the results of previous works^iS when only the spin- flip 
scattering is considered. 

We should point out that our formalism applies also 
to metals, by considering the case ^ ^ M. In this case, 

the last term of Eq. can be neglected and the effect 
of the spin-conserving scattering through t^* becomes ir- 
relevant. 

When the pure spin current is included, we found ad- 
ditional contribution due to the interplay of the spin cur- 
rent and the SOC, since we have 



n 



k X 



(0) 



(23) 



with the spin current associated term s^^ defined ac- 
cordingly. The origin of s^^ can be understood as fol- 
lows. The anisotropic spin polarization sf^y^ arising from 
the pure spin current rotates around the SOC effective 
magnetic field hk, which is also anisotropic. This pre- 



cession finally results in an isotropic spin polarization 
s^'^ = Tr J2u '^k X sfl. in the presence of spin relaxation. 
This term contributes to the spin polarization of the itin- 
erant electrons along the direction of the magnetization, 
i.e., sl^, thereby modifies the Gilbert damping term by 

The additional Gilbert damping due to the spin cur- 
rent found here is different from the enhancement of the 
damping in the spin pumping systems, where the exis- 
tence of the interface is essential.* In other words, what 
contributes there is the divergence of the spin current, as 
is understood from the continuity equation for the spin, 
indicating that the spin damping is equal to V • js -I- s (s 
is the total spin density). In contrast, the damping found 
in the present paper arises even when the spin current is 
uniform if the spin-orbit interaction is there. 

In summary, we have shown that the spin-conserving 
scatterings in ferromagnetic semiconductors, such as 
the electron-electron, electron-phonon and electron-non- 
magnetic impurity scatterings, contribute to the Gilbert 
damping in the presence of the SOC because of the in- 
homogeneous broadening effect. We also predict that a 
Gilbert torque arises from a pure spin current when cou- 
pled to the spin-orbit interaction. 
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